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Abstract 

A nonlinear frequency response based adaptive vibration controller is proposed for a class 
of nonlinear mechanical systems. In order to obtain the nonlinear Frequency Response Function 
(FRF), the convergence properties of the system are studied by using the convergence (contrac¬ 
tion) theory. If the system under consideration is: 1) convergent, it directly enables to derive 
a nonlinear FRF for a band of excitation inputs, 2) non-convergent, first a controller is used to 
obtain the convergence and then the corresponding FRF for a band of excitation inputs is derived. 
Now the gains of the proposed adaptive controller are tuned such that a desired closed-loop fre¬ 
quency response, in the presence of excitation inputs is achieved. Finally, a building structure 
with nonlinear cubic stiffness and a satellite system are considered to illustrate the theoretical 
results. 
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1. Introduction 

Mechanical vibrations are present in countless real-life situations, where the mechanical sys¬ 
tems exhibit oscillations when subjected to certain excitations. Most often these vibration phe¬ 
nomena are highly undesirable, which may even cause damage to the system itself. Vibration 
analysis and control are an active, vast, and growing research area, due to its practical impor¬ 
tance and issues that arise in both linear and nonlinear system designs. The vibration control is 
concerned with the prediction and controlling of these undesired oscillations. The vibration can 
generally be controlled by adding controlling devices like dampers, isolators, and actuators to 
the system. These devices are added in such a way that the system’s properties are modified to a 
desired one (^j. Most of the vibration control methods work based on time-domain techniques. 
These methods lack in describing how a closed-loop system respond to the input excitation at 
different frequencies and magnitudes. Since the vibration is characterized by its frequency (or 
frequencies), amplitude, and phase, it is important to study the frequency response of these sys¬ 
tems. 

Frequency domain techniques for linear systems have led to significant progress in analysis, 
modeling, and controller design Mi. In reality, mechanical systems posses many physical 
properties such as material property, geometric nonlinearity, damping dissipation, and even due 
to boundary conditions, which lead to nonlinear vibration problems ||3l- For that reason, these 
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nonlinear systems do not posses simple oscillations as defined in the case of linear systems. 
These nonlinearities can result in some complex phenomenon, like jumping, chaos, secondary 
resonance, bifurcation, etc. 0. In those instances, the classic linear frequency response analysis 
tools are insufficient to describe the behavior of the nonlinear system adequately. 

The methods such as describing function (DF) lack in accuracy due to its approximation 
scheme ||3l. Nonlinear FRF such as the Generalized Frequency Response Function (GFRF) 0 
is limited to second order due to its multi-dimensional characteristics. An extension of the GFRF, 
termed as Output Frequency Response Function (OFRF) was proposed in which repre¬ 

sents the relation between the system parameters and frequency response using finite Volterra 
series. Flowever, this technique fails to detect some of the nonlinear phenomena such as the sub¬ 
harmonics, jumping, etc. Ml ill . In lll2ll . it has been shown that an FRF can be found for a class 
of nonlinear systems termed as convergent systems. The convergence property implies that the 
system trajectories tend towards the unique bounded solution. If the system under consideration 
is convergent, it signifies; 1) the system is stable El, 2) an FRF can be obtained due to the 
existence of a unique steady-state solution El. Compared to the DF, GFRF, and OFRF, this 
FRF gives an exact frequency response via numerical or experimental approach. 

The main objective of this paper is to introduce the potential of convergence analysis and 
nonlinear FRF in vibration control problems. This paper considers a nonlinear FRF based vibra¬ 
tion analysis and control of a class of nonlinear mechanical systems. In order to derive the FRF, 
hrst the convergence properties of the system have been established theoretically. In the case of 
non-convergent systems, a feedback control is designed to achieve convergence. In terms of con¬ 
trol, the controller gains are adapted based on the FRF of the system derived within the excitation 
band of interest, which assures a satisfactory performance over that band. The control scheme 
developed here is applied to a building structure with cubic stiffness nonlinearity and to a satellite 
system, and the corresponding dynamic responses of both the controlled and uncontrolled cases 
are numerically evaluated. 


2. FRF of Nonlinear Convergent Systems 

For a system to be satisfactory, it is necessary to analyze its stability. In general a system’s 
stability is analyzed by examining, whether the equilibrium points so determined are stable. The 
convergence analysis, inspired by fluid mechanics, is the extension of the stability properties 
of asymptotically stable linear time-invariant systems. Unlike the Lyapunov stability theorem 
which defines the stability with respect to the equilibrium points, the convergence in a convergent 
system implies that the state trajectories with different initial conditions will converge to a unique 
bounded solution El- 


Definition 2.1. Let a dynamical nonlinear system be described by the differential equation 

X = fix, t) 


( 1 ) 


with X e R" is the state vector, t e R+, and / : M" x R+ —> R" is a smooth nonlinear function. 
For the above system, a region of convergence (or contraction region), X is defined where the 
system's Jacobian matrix, J(x) — df jdx is uniformly negative definite. 

The convergence properties of the system ([TJ can be verified by performing a coordinate 
transformation on 7(x). The resulting generalized Jacobian is dehned as 




-H Y7(x))t- 
2 


1 


( 2 ) 



where T (x, t) is a uniformly invertible square matrix. If Sf is uniformly negative definite 

J < -dtnaxf, 'ixeX d R", Vf 6 K (3) 

where dmax is the largest eigenvalue of Sf, then the transformed system (|2]l is convergent, which 
implies that all the solutions of the original system ([TJ converge exponentially to a single trajec¬ 
tory, independently of the initial conditions. If S' is negative semi-definite, then the system is 
semi-convergent under some mild conditions similar in Barbalat’s lemma, this implies that the 
solutions converge each other asymptotically. The global convergence or semi-convergence is 
achieved when A' = R". 

FRF is the characteristics of a system that describes its response to an input as the function 
of frequency. Consider a nonlinear time-invariant system, which is forced by the input excitation 
W £ 'W 

y = gix) ^ ’ 

where y is the system output. 

Definition 2.2. If the system m is convergent, then there exists an uniformly bounded steady- 
state (UBSS) solution, if for any p > 0 there exists cr > 0 such that for any input w £ TV the 
following implication holds: 

\w\< p Vf e R \x„\ < cr Vf 6 R (5) 

where x„, is the steady-state solution, which depends on w. For the convergent systems, this 
bounded solution x„ is unique and lim,^oo \\xwit) — x„(f)ll = 0/or any xq £ X and hence the 
system (0 is exponentially stable /fH/- 

Definition 2.3. If the system 0 is convergent with UBSS property fora certain class of harmonic 
inputs w(t) — a sin(a»f) £ TV, then there exists a nonlinear function a : R^ —> R" such that 

x„(t) ;= a(vi, V 2 ,w) (6) 

where Vy — a sin(a»f) and V 2 — a cos(iot). The nonlinear function a (vi, V 2 , co) is known as the state 
frequency response function and the function g (a (vi,V 2 , co)) is known as the output frequency 
response function, which relates the harmonic input to the corresponding steady-state output 


The output of the forced system at steady-state can be expressed as y^„(t) - g(x„(t)), which 
will have the same period of the input signal w, but not necessarily sinusoidal. Now the output 
response of the system for various amplitude and frequency inputs can be represented using an 
amplification gain ja (ai), which is the ratio between the maximal absolute value at steady-state 
and the corresponding input signal amplitude, so that 


1 

ya,(jj — 

a 


sup |g(a(vi,V2, w))| 


|y>v(0| 

a 


(7) 


Analyzing the above function, one can find the critical input amplitudes and frequencies for 
the system and an appropriate controller can be designed to bypass any undesirable effects such 
as the vibration. 
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3. Vibration Control of Nonlinear Mechanical Systems 

There are three basic elements in a mechanical system: mass which stores kinetic energy, 
spring which stores the potential energy and the damper which dissipates the energy. When 
an external force is applied to a structure, it produces change in its displacement, velocity, and 
acceleration. 

Let us consider a nonlinear multi-degree-of-freedom mechanical system of the following 
form 

Mq + Cq + Kq + ^(q, q) — Aw (8) 

where M, C, and K e Q are the positive definite matrices correspond to the mass, damping, 
and stiffness respectively, q, q, and q e R."'' are the relative acceleration, velocity, and displace¬ 
ment vectors respectively, A 6 R”? denotes the influence of the input excitation force w on the 
system, and ^{q, q) e R"'' is the system nonlinearity vector. 

3.1. Problem Formulation 

In vibration control applications, the system dynamics were modified favorably by adding 
passive or active devices. Most of the passive devices can be tuned only to a particular structural 
frequency and damping characteristics. But the nonlinearities in the systems cause variations in 
its natural frequencies and mode shapes. Since the passive devices cannot adapt to these struc¬ 
tural response changes, it cannot assure a successful vibration suppression. These shortcomings 
can be overcome by using an active control system, where the system’s output response is mea¬ 
sured using sensors and an appropriate control force, calculated by a pre-assigned controller is 
used to drive the actuators for suppressing the unwanted structural vibration. The mechanical 
system with an active control system can be represented as, 

Mq + Cq + Kq + ^(q, q) = Aw + Fu (9) 

where u e R"" is the control signal generated by a control algorithm by means of which the 
system can be controlled and T e is the actuator location matrix. 

Now let us consider a feedback controller of form 

M = -&x (10) 

where © e R"“^" is the state feedback controller gain matrix and x — ^q^ q^^ . The control 
objective is to find a frequency response based adaptive control rule for the controller gain 

© = T'(D (11) 

where 'K is the FRF of the system, such that a desired vibration attenuation by means of a 
specified frequency response is achieved for the closed-loop system (|9]l. 

3.2. Vibration Control of Convergent Systems 

In this section, first we will establish the convergence property of a nonlinear mechanical 
system ([8]l, which will help us to derive the FRF of the open-loop system. Next we will derive the 
range of controller gains for which the closed-loop system is convergent, hence 0 e Q© c R®, 
where Q© is a compact set. Finally, based on the system’s FRF, an adaptive algorithm of form 
(fTTll will be presented to tune the controller gains 0 within a convergence region such that the 
vibration is minimized. 


*In the case of second-order mechanical systems n = 2n^. 
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3.2.1. Convergence Analysis 

Let us consider the multi-degree-of-freedom mechanical system shown in ([8]l, with an odd 
polynomial nonlinearity of form 




( 12 ) 


where 


p 

^ bigp+iq^‘’*\ bigp+i > 0, i = 1, n^ 

p=i 


with p as the highest order of the nonlinearity. By letting q - xi e and q - X 2 & K"’, the 
corresponding second order model represented can be written as a set of first order differential 
equations 


Xi 


X2 

X2 


[Cx 2 + Kxi + O/xi) - Aw] 


(13) 


which is equivalent to the nonlinear model (|4|, where the nonlinear vector functions f{x, w) are 


/l(Xi,X 2 ,w) 


■^2 

/2(X1,X2,W) 


[Cx 2 + Kxi + <l>(xi) - Aw] 


(14) 


Theorem 1. The class of mechanical system with odd polynomial nonlinearity, represented in 
is convergent. 

Proof. The system (0 can be represented in terms of virtual displacement 6x as 10/ 


6x - J(x)6x 


(15) 


The Jacobian matrix of the system represented in m is 


J{x) = 


On^xn, 

-M-l [/T + Oy/^i)] 


-M-^C 


e R 


nxn 


(16) 


where <I>y(xi) = <I)y(xi) = ^ 0- Generally, it is difficult to obtain a uniformly negative 

definite Jacobian matrix directly for the system represented in the form m- However, a trans¬ 
formation can always be done to prove the negative definiteness of the Jacobian matrix. Consider 
the following coordinate transformation 


6z = Tdx 


where the transformation matrix T is 


T = 


bn,iXn,i 0«,xn, 
bn,,xng ^n,xn. 


then 


dzi 


6x\ 

6Z2 


(5xi + 5x2 


(17) 


(18) 


(19) 
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Using m the generalized Jacobian JJ can be calculated as 


J 


rj(x)r-^ 

-In,xn, - M-^ [K + Oj(zi) - C] 


n -M-^C 


( 20 ) 


Now the transformed system is 


6z - fjdz 


( 21 ) 


The squared distance between the trajectories of the virtual system can be calculated as 


vll&ll = 26z^6z 
at 

= 26z^J6z 

< ( 22 ) 


In order to achieve exponential convergence, ff must be negative definite, which is estab¬ 
lished since ffu < 0 and K + <l)y(zi) > 0 (/T > 0 and <l*y(zi) > 0). Hence, the solutions of the 
virtual system ( 1271 ) converge exponentially to zero with a rate (Jj, which implies that the 
solutions of the actual system (H also converge to each other. 

3.2.2. PD Controller Design 

PD control may be the simplest controller for the vibration control system, which provides 
high robustness with respect to uncertainties. It has the following form 

u — -&p e - &de (23) 

where &p, 0y e are symmetric matrices correspond to the proportional and derivative 

gains, respectively, e - i^q- q‘^^ corresponds to the position error, e - i^q- corresponds to 
the velocity error, and q‘^ is the desired position. In active vibration control case, the references 
are q^ - ijf = 0, hence (|2^ becomes 


u = -®pq - &dq (24) 

The closed-loop system (|9]l with the PD controller (l24l i is 

Mq + Cq + Kq + O)^') = Aw - Opq - &dq (25) 

which can be written in the state-space form as 


ri 

X2 


X2 


-M ^ [ (C H- @d)x 2 -H (tT -H ©p)ri -I- <t(xi) - Aw ] 


(26) 


The range of gains for which the system is convergent is defined in the following theorem. 


^In the case of full-state feedback control In that case F is an Identity matrix of dimension itg, hence ignored 

in the later part of this paper. 
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Theorem 2. Consider the class of nonlinear system m with an external excitation w controlled 
using the control law m- If we choose the control gains such that 0 = ©^j 6 Q© c K®, 

then the state trajectories of the closed-loop system & lies within the region of convergence 
"xw & X and hence exponentially stable. 

Proof. The Jacobian matrix of the closed-loop system ( 1261) is 


J{x) = 


On,xn, ^n,XK, 

-M-^Kk -M-^Kc 


(27) 


where Kk = /T + ©^ + <]> y(xi) and Kc — C + @d- The generalized Jacobian ff can be now found 
as 


J 


rj(x)r-^ 

~In^xn^ 

-In,xn, - M-i [Kk - Kc] 


-M-^Kc 


(28) 


Since Kk > 0, the matrix ff is negative definite. From the above discussions and from ([5]), it can 
be concluded that for a bounded input excitation |w| < p, the controller with positive proportional 
and derivative gains will result an exponentially stable solution, hence 'x„ e X. 


3.2.3. Nonlinear Frequency Response based Adaptive Controller Design 

Adaptive techniques are widely used for active vibration control applications. Traditional 
time-domain based adaptive schemes update the controller gains based on the error at that mo¬ 
ment. However, they cannot assure a satisfactory vibration attenuation for a given band of exci¬ 
tation. Here, designing of an adaptive control algorithm based on the system’s FRF is discussed. 

The amplification gain ja (w) of the system, for a range of magnitudes (a e (a, a) e K'') and 
frequencies (a> e ZZi) e M^) can be represented as 


To 




Ya,^ • ■ • Ya,to 


(29) 


The above matrix can be considered as the open-loop FRF matrix of the system. Now Fq 
can be analyzed in order to get a knowledge about the critical magnitudes and frequencies of the 
excitation input, at which the system possess a larger amplification gain. One way to evaluate 
the FRF matrix is by finding its Frobeinus norm (F-norm), which are sensitive towards its each 
elements. F-norm of a matrix A e can be calculated as 


IIAIIf - 


A 




(=1 y=l 




>0 


(30) 


For the critical input excitations, the F-norm of the open-loop system FRF \\To\\p will be 
higher. A PD controller can be used to reduce the system peaks at those critical points. The 
proportional and derivative gains provide virtual stiffness and damping to the closed-loop sys¬ 
tem, respectively. The damping term helps in reducing the resonance peaks by dissipating the 
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vibration energy. Once we increase the A^ax of the closed-loop system becomes more nega¬ 
tive and as a result the convergence rate (e'*™") increases accordingly. On the other hand @p can 
be tuned such that the system resonance frequency can be shifted beyond the frequency band of 
potential excitation. The stiffness value must be chosen carefully without causing any undesir¬ 
able effects at other excitation inputs. This phenomenon is termed as waterbed effect, where the 
reduction of the magnitude response of the closed-loop system in one frequency range will result 
in the increase of the magnitude response in some other frequency range Q.' By analyzing the 
system’s closed-loop FRF, the controller gains can be tuned in order to achieve maximal vibra¬ 
tion attenuation at the critical excitations without causing undesirable effects at other operational 
frequencies. 

The steady-state output of the closed-loop system for a particular value of proportional and 
derivative gains can be represented as y„(t, 0). The amplification gain of the closed-loop system, 
denoted by ■ya.oj (©), satishes the following relation 

|yw(f,0)| = r«,-(0)N OD 

Flence the peak vibration output of the system can be attenuated by minimizing the amplifi¬ 
cation gains. Using the amplihcation gains obtained under a range of excitation, the closed-loop 
FRF matrix, can be constructed similar to (1291) . The control objective is to minimize the FRF 
magnitude such that 


\\r@\y = 


\l/2 


y! y (0 )/j 


V/=l j=l 


<s< wnwF 


(32) 


where 6 is the acceptable vibration range. For the ideal vibration attenuation case 6 may be 
chosen as zero, which means 'ya,oj(&) = 0, due to the property HlFell/r = 0 'Fe - 0. 

However, in practice it is not possible to remove the entire vibration from the system due to the 
waterbed effect and limitations of the control devices. Our goal is to choose a d according to the 
practical situations and to achieve it by using an adaptive scheme. 

Here, the design requirement is that each increment in the controller gains at least minimizes 
the F-norm of the FRF matrix. The controller gains are adapted by using the following algorithm 

A0,-40,.^i- 0,. = Feller®,,11^ (33) 


where 

Te-F^ 

is the adaptation step size, 

WreA^l, 0 
0 WreAMl: 

and Si = (Ii)||^ - , \\r@,i (T 2 )||^ - dxjj is the error between the F-norm of the closed- 

loop FRF at /-th iteration and the desired range 6. If the gains are chosen such that 



To,, 0 
0 r. 


> 0 


0, if 0, > 0niin 

©tnin if 0/ < ©min 
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( 34 ) 








where ©min = [©pmin ©i/min] IS used to assure that the controller gain 0; is always positive, then 
the solutions of the closed-loop system stays within the convergence region. 

Equation (1331) is used to calculate the new controller gains, according to its previous FRF 
matrix. When the error e is positive, the gains will increase and for negative e the gain decreases. 
Based on the e, the gains will adapt over each iteration until a satisfactory vibration attenuation 
is achieved for a band of excitation, that is HlFell/r = d, hence s, ^ 0 as i —» oo. Moreover, the 
projection operator (l3^ will assure that the closed-loop system (l25T l with the adaptation scheme 
(l33T l operates within the convergence region. 

3.3. Vibration Control of Non-Convergent Systems 

In the previous case the open-loop configuration of the system under consideration is al¬ 
ready convergent. But in practice, there exist mechanical systems which are not exponentially 
convergent. Let us consider a nonlinear mechanical system of form 

Hiq)q + Ciq,q)q = T (35) 

where H{q) e is the symmetric positive definite inertia matrix, C{q, q)q e K"« is a non¬ 

linear vector function of q and q due to the Coriolis and centripetal effects, and t e R"’ is the 
actuator input torque. Since, for the above system the gravitational force term {dUjdqf corre¬ 
sponds to the potential energy U is absent, and the corresponding total energy of the system E, 
is equal to its kinetic energy, which is 

E = ^fHq (36) 

and the corresponding variation of the energy is 

E - q^Hij + -q^Hq 

= i-Cq + t) + ^f^Hq 

= ^f[H-2Cq) + fT (37) 

Considering t = 0 and by using the skew-symmetric property of the system dLSl) such that 

[H(q) - 2C{q, q)]q ^ 0 (38) 

yields E - Q, which implies the total system’s energy is constant due to the conservation of 
energy. This is due to the fact that the matrix C{q, q) does not contain any dissipative forces. 
Now the non-conservative force, like the input torque (E - q^r), might be designed such that 
system trajectories converges exponentially. In this section, we discuss a method to provide 
exponential convergence by means of a feedback control loop. 

3.3.1. Obtaining Convergence 

One way to achieve the system convergence is by introducing a virtual stiffness and damping 
term by means of a standard PD controller. Applying the PD control law in (1351) yields 


H{q)q + C{q, q)q + Koq - -Kpq 
9 


( 39 ) 


Now let us define a local coordinate system z defined as 

i^H6v + C6q + Kodq = -Kpq (40) 

The virtual displacement of the above system can be written as 

dz = HSv + C5q + KD6q (41) 

and the corresponding dynamic variation is 

6z = -Kp6q (42) 


Now the convergence can be established as follows, 


^ [dv^Kp^dv + 6q^H6q\ 


-26v^5q + 25q^ {H6v + C6q) 


-26q^KD6q<0 


(43) 


Using Barbalat’s lemma it can be shown that both the dv and 6q will converge to each other 
asymptotically, which implies the semi-convergence of the original system (ITST i. For that reason 
the system FRF cannot be obtained, which requires that the system states reach steady-state 
exponentially. 

Now let us consider an energy-based controller of form OlSll 

T = H{q)q'' + C{q, q)q'' - K^iq - (f) (44) 

where is a symmetric positive definite gain matrix and the reference velocity 

q'' - if - KyC, Ar > 0 
and the corresponding tracking error r is 


r - e + ArC - q- f (45) 

The virtual system is 

T - H(q)z +C(q,q)z - Kriq - z) (46) 

To examine the convergence properties, the virtual distance can be calculated as 

^ [dz^Hdz] = -26z^ {Hq - Hif) + dz^Hdz 

= -2Sz^ i-Cq + T- Hq^) + 6z^H6z 

= -26z^ (-Hq'' - Of t) -H [h - 2C) 6z (47) 

Once again using the skew-symmetric property (l38T l of the system dTSl l and the controller 
(l44li . yields 

^ [dz^Hdz] = -25z^Kr6z < 0 (48) 

which indicates that q converges to if exponentially. By noting that the trajectories of dTSl l and 
(l44l l are particular solutions of the virtual system (l46l l its convergence signifies that q tends to q^ 
exponentially. Now an active vibration controller can be used to minimize the vibration caused 
by excitation inputs. 
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3.3.2. PD Controller Design 

In order to analyze the convergence property the PD controller law can be rewritten as 


u - -Orf - -&r (e + A^e) (49) 

where = ©rA^, &d - ©,. Now the energy-based controller (l44l) along with PD controller 
(l49l l can be represented as 

T = H(q)q'- + C(q, M " K.iq - q’’) - &Aq - ql (50) 

Theorem 3. The nonlinear mechanical system 051 ) with energy-based controller is conver¬ 
gent for the PD controller l \49h such that © 6 Q© c R®. 

Proof. The virtual system is 

T = H(q)z -t- C(q, q)z - KAq - z) - ©Aq - z) (51) 

Using the controller m the virtual distance can be calculated as 

^ [dz^Hdz] = -2&^ (Kr + &r) & (52) 

which signifies that q converges to q' exponentially, therefore, exponential convergence of q to 
qd is guaranteed. 

Since the controller in (l49l l has the structure and gain ranges similar to that of (l24l i. we are 
in a position to implement the frequency response based adaptive controller proposed in Section 
l3.2.3l to the system dlST l with controller (l44l) . 

4. Results and Discussions 

The frequency response based vibration control algorithm proposed in this paper can be used 
in a variety of vibration control applications. In this section, we have considered two important 
vibration control applications: 1) a building structure with cubic stiffness, which is a convergent 
system, 2) a satellite system, which is a non-convergent system. The performance of the proposed 
algorithm was evaluated via numerical simulations. 

4.1. Application to Building Structures 

Protection of large civil structures and human occupants from natural hazards like an earth¬ 
quake and wind is very important vibration control application. Many attempts have been made 
to introduce advanced controllers for the active vibration control of building structures 111]. Most 
of these methods work based on time-domain techniques. In practice, the building structures 
behave nonlinearly under large deformations, which can happen during strong seismic events. 
During these situations the well established linear frequency response analysis tools cannot be 
applied. 

The cubic nonlinearity have been chosen here, which have received considerable interest in 
literature. In order to present the main idea, let us consider a single-degree-of-freedom mechan¬ 
ical system with a cubic stiffness element ike) 

mq cq kq kcq^ — w 
11 


( 53 ) 




where the parameters are set as: m — 1kg, c = 0.4Ns/m, and k - kc - 36N/m. The control 
programs were operated in Windows 7 with Matlab 8.0/Simulink. All the control actions were 
employed at a sampling period of 10ms. The magnitudes and frequencies of the input excitation 
are a(N) e (0.5,6) and w(rad/s) 6 (3,9). The block diagram of the implementation of the 
proposed algorithm is shown in Figure[T] 



Figure 1: Block diagram of the proposed active vibration control system. 


The convergence of solutions for different initial conditions are shown in Figure From this 
plot, it can be seen that the solutions for different initial conditions (ji:i(0) = -3 ,jc2(0) = 3,Xi(Q) = 5) 
converge to the equilibrium point. Due to this convergence property, the FRF from the steady- 
state response can be derived without considering different initial conditions. 



Figure 2: Convergence of trajectories for different initial conditions. 


Figure [3] shows counter plots of the open-loop FRF of the mechanical system without and 
with cubic stiffness, respectively. For the linear case the resonance frequency is = 6rad/s. The 
linear case shows that the magnitude of the external excitation has no effect on the amplification 
gain, hence the resonance frequency of the system. But in the nonlinear case, the amplification 
gain is a function of both the magnitude and frequency of the input excitation, which causes a 
shifting in its resonance frequency for different values of a and o). The active control scheme can 
adapt to these nonlinear effects. The closed-loop system is 


mij + cq + kq + k^q^ — w — 9pq — Odq (54) 

The effects of adaptive PD controller on the performance of the active vibration control sys¬ 
tem are investigated. The adaptive algorithm parameters are set as: Fg^ = 120,Fg^ = - 

0.5,- 3, and ^min = [0.001 0.001]. The controller gains was started from a minimal value 
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(a) Linear case {kc = 0) 
Figure 3; Open-loop FRF 
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(b) Nonlinear case {k^ = 100) 
of the mechanical system. 



(a) Evolution of adaptive gains 
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(b) Evolution of adaptive error 


Figure 4; Evolution of adaptation parameters over each iteration. 


0ta - Omin, and allowed to adapt using (l33t after each set of FRF is calculated. Figure|4l(a) shows 
the controller gain adaptation process of the closed-loop system (l54l i. The corresponding error 
function s is shown in Figure|4l(b). Figure|5]compares the FRF of the system without and with 
control, respectively. The F-norm of the system for position measurement, without (||!To(^i)IIf) 
and with control (||'Ts(Ti)||/7) are 1.42 and 0.5, respectively. The controller response is obtained 
using the final values obtained from the adaptive scheme, i.e. Op = 7.1N/m, 0^/ = 2.6Ns/m. 
Figure |6l(a) shows the time response of the displacements for both controlled and uncontrolled 
cases and the corresponding PD controller output signal is shown in Figure |2(b). 

From these simulation results, it is clear that the proposed algorithm assures an effective vi¬ 
bration attenuation over the band of excitation. Nonlinearities such as cubic stiffness can cause 
subharmonics in the system response. Since the F-norm of the system FRF is sensitive to any 
resonance in the excitation band, the iteration process is repeated until a stopping criterion is 
satisfied. Consequently, the new algorithm can effectively enhance the control performance over 
the given range of inputs. Furthermore, the stability analysis is much simpler due to the conver¬ 
gence analysis, which is obtained by performing simple coordinate transformation. Note that the 
controller updates its gains based on the measured input-output (amplification gain) relation of 
the system, which means that the exact knowledge of the system parameters is not required for 
the gain adaptation. 


4.2. Application to Satellite Attitude Control 

Momentum type actuators such as Reaction Wheels (RW) are widely used as attitude control 
actuators in spacecraft for orbital maneuvering. Unfortunately, these actuators are one of the 
main source of on-board vibration lllhll . which is caused by the static and dynamic imbalances 
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Figure 5: Uncontrolled and controlled {Op - 7.1,0^ = 2.6) FRF of the system. 




(a) Uncontrolled and controlled displace-(b) Control signal generated by the PD 
ments of the system. controller. 

Figure 6: Time response of the closed-loop system. 


in the RW assembly. This can be critical for the high precision space applications such as high- 
sensitivity imaging and astrometry. The dynamics and kinematics of a satellite system can be 
modeled as 

Hus + S {UsjHus - T (55) 

q - JstOs (56) 

where H e is the inertia matrix, a>s e is the angular velocity, S {tOs)H is the angular 
momentum with S (.), a skew-symmetric matrix representing the vector cross product, t 6 
is the torque applied to the satellite system, g' e R^ is the satellite attitude vector, and 6 
M3x3 

: tUj —> 4' is the Jacobian matrix, all expressed in the satellite body frame. In this paper, the 
Modified Rodrigues parameters were used to represent the kinematic equations of the satellite 
system. 

By considering q and q as the state-space coordinates and using (l56]) . the equation of motion 
of the satellite system (l55]) can be written in the Lagrangian form as 

Hs{q)'q -I- Csiq, q)q = T, (57) 

where 

HM) = 

CAq,q) = jfHJ-^Jj-^-jfS{cOs)HJ-^ 
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Since the system (l57l i verifies the structure and skew-symmetric property ifP^ of dlSl l. it 
satisfies the theoretical analysis presented in the Section [53] 

The rotational elements of the RW generates periodic disturbances, which can be modeled as 
a series of discrete harmonics Q 

h 

Wrw - sin (2jThiD.t + ai) (58) 

i=l 

where h is the number of harmonics, A, is the amplitude coefficient of the /-th harmonic, Q is 
the wheel speed, is the /-th harmonic number and a,- is a random phase. Since our disturbance 
model w is similar to w™ by choosing a - A,Q^ and oj = 2nhiQ., the proposed method can 
be a potential candidate in designing a vibration controller for the spacecraft systems with RW 
assembly. Now the satellite model with the disturbance Wrw can be represented as 

Hus + S (a>s)Ha>s - t + u + w™ (59) 

where u is generated by the proposed FRF-based adaptive control algorithm for minimizing the 
vibration signals caused by the RW assembly. The block diagram of the implementation of the 
proposed algorithm is shown in Figure [7] 


Controller 



Figure 7; Block diagram of the active vibration control system for satellites. 


Combining equations (l50l l and (fSTT i yields, 

Hs(q)q + Cs(q, q)q = HiqYq'' + C(q, q)q’' - (Kr -I- &r)(q - q') (60) 

where q'' = q‘^ - Are, and the corresponding closed-loop dynamics can be written in the following 
simple form 

Hs(q)r + (Cs(q,q) + Kr+ &r)r-0 (61) 

In order to implement the controller (l44li . we need to have the information regarding the 
parameters of the satellite. The satellite parameters utilized in ifisll are used here to illustrate 
the numerical results. The initial conditions are set as q^ - q^ - [0 0 0]^. The satellite sys¬ 
tem is excited using sinusoidal disturbances, and the corresponding position and velocity error 
((F© (e) and 'Tq at steady-state were used to obtain the FRF. The magnitudes and frequen¬ 
cies of the disturbance w™ for different RW velocities can be estimated as proposed in M- 
During the tuning stage, the desired input q^ is selected such that q and q reaches steady- 
state. In this case, the satellite desired attitudes are chosen to be q‘^ - [1 0.5 0]^. The adap¬ 
tive algorithm parameters are set as: Fe^ = diag[10 10 10], F©^ = diag[1000 1000 1000], 6j - 
[0.2 0.1 0.1], 6^ = [0.02 0.02 0.02], and =0.001. 
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Figure 8 : Evolution of adaptive gains. 



(a) Proportional gain error (b) Derivative gain error 

Figure 9; Evolution of adaptive errors. 


The adaptation of the proportional and derivative gains based on the FRF are shown in Fig- 
ureOand the corresponding evolution of the errors, (||lF 0 ,, (F)||^ - d?) and (^)|| 

shown in Figure|2(a) and (b), respectively. It is seen that, by using the controller (l44li the position 
and velocity error converges to zero quickly. However, the controller performance degrades in the 
presence of any excitation inputs. Figure[TO]shows the positional and velocity error of the satellite 
in the presence of excitation input, without using the adaptive controller (u — 0). From these dia¬ 
grams we can see that the controller requires some time interval to bring the error to zero. Figure 
QTIshows the positional and velocity error of the satellite in the presence of excitation input, using 
the FRF based adaptive controller (@p = diag[0.93 1.98 1.53], and = diag[10.60 5.24 5.13]). 
From these diagrams, it is seen that the FRF based adaptive controller provides a suitable pro¬ 
portional and derivative gain such that the closed-loop system has less sensitivity towards the 
excitation force, hence induces the error caused by the excitation to converge quickly and as a 
result a better attitude control performance is achieved. 

5. Conclusions 

This paper has described an adaptive vibration control scheme for a class of nonlinear me¬ 
chanical systems, subject to external excitation by using the nonlinear FRF technique. The non¬ 
linear system possesses complex behavior for different excitation inputs. A simple and effective 
tool for vibration analysis and controller design has been proposed. Theoretically established 
that a mechanical system with an odd polynomial nonlinearity is convergent and hence a FRF 
can be derived. Furthermore, the stability is assured due to its convergence property. For the 
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(a) Position error (b) Velocity error 

Figure 10: Without FRF based adaptive controller. 



(a) Position error (b) Velocity error 

Figure 11: With FRF based adaptive controller. 


systems that are not convergent, a controller is used to establish the convergence. The main ad¬ 
vantage of the frequency-domain approach presented in this article is that a satisfactory vibration 
attenuation for a band of excitation is assured. Finally, the proposed scheme is applied to the 
active vibration control problem of building structures with cubic nonlinearity and satellite sys¬ 
tems. From the numerical studies, it is observed that the procedure is successful in suppressing 
the vibration. 
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